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Abstract 

Nuclear pairing properties are studied within an approach that includes the quasiparticle-number 
fluctuation (QNF) and coupling to the quasiparticle-pair vibrations at finite temperature and an- 
gular momentum. The formalism is developed to describe non-collective rotations about the sym- 
metry axis. The numerical calculations are performed within a doubly- folded equidistant multilevel 
model as well as several realistic nuclei. The results obtained for the pairing gap, total energy and 
heat capacity show that the QNF smoothes out the sharp SN phase transition and leads to the 
appearance of a thermally assisted pairing gap in rotating nuclei at finite temperature. The cor- 
rections due to the dynamic coupling to SCQRPA vibrations and particle-number projection are 
analyzed. The effect of backbending of the momentum of inertia as a function of squared angular 
velocity is also discussed. 
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I. INTRODUCTION 



Thermal effect on pairing correlations has been extensively studied within the Bardeen- 
Cooper-Schrieffer (BCS) theory fl] at finite temperature T (FTBCS theory). The FTBCS 
theory predicts a destruction of pairing correlation at a critical temperature ^ 0.568A(0) 
[A(0) is the pairing gap at zero temperature], resulting in a sharp transition from the super- 
fluid phase to normal one (the SN phase transition) in good agreement with the experimental 
findings in macroscopic systems such as metallic superconductors. However, the BCS theory 
is valid only when the assumption on the quasiparticle mean field is good, i.e. when the 
difference between the pair correlator P^^ = al.^al_^ and its expectation value (pI^) is small 
so that the quadratic term {pI^ — (-Pfco-))^ negligible, where a|,^ is the operator that creates 
a particle with angular momentum k and spin a. For small systems such as underdoped 
cuprates, where the coherence lengths (the Cooper-pair sizes) are very short, the fluctuations 
{pI^ — {Pla)y longer small, which invalidate the quasiparticle mean-field assumption, 

and break down the BCS theory. As the result, the gap evolves continuously across Tc, and 
persists well above Tc Q]. 

Various theoretical studies have been undertaken in the last three decades to study the 
effects of thermal fluctuations on pairing in atomic nuclei. Pioneer papers by Moretto [Sj] 
employed the macroscopic Landau theory of phase transition to treat thermal fluctuations 
in the pairing field as those occurring around the most probable value of the pairing gap. 
The results of calculations within the uniform model carried out in Ref. 3j] show that 
the average pairing gap does not collapse as predicted by the FTBCS theory, but decreases 
monotonously with increasing T, smearing out the sharp SN phase transition. This approach 
was later used by Goodman to include the effects of thermal fluctuations in the Hartree-Fock- 

n 

Bogoliubov (HFB) theory at finite temperature The static-path approximation (SPA), 
which takes into account thermal fluctuations by averaging over all static paths around 
the mean field, also shows the non-collapsing pairing gap at finite temperature 5|, l6| . The 
recent microscopic approach called the modified BCS (MBCS) theory is based on 

the secondary Bogoliubov transformation from quasiparticles to the modified ones to restore 
the unitary relation for the generalized particle-density matrix at T 7^ 0. The MBCS theory, 
for the first time, points out the quasiparticle-number fluctuation (QNF) as the microscopic 
origin that causes the non-collapsing thermal pairing gap in finite small systems. The 
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predictions of these approaches are in quahtative agreements with the experimental findings 



ot paMng gaps and heatcapacMes measured .n unde.doped cup.ate. H and extracted fro,n 
nuclear level densities [10] • 

While quasiparticles are regarded as independent in all above-mentioned approaches, 
the recently proposed FTBCSl theory with corrections coming from the QNF and self- 
consistent quasiparticle random-phase approximation (SCQRPA) at finite temperature 11 1 
calculates the quasiparticle occupation numbers from a set of FTBCSH-SCQRPA equa- 
tions. Within this approach, which is called the FTBCSH-SCQRPA and is an extension 
of the BCSH-SCQRPA developed in Ref. [isl to finite temperature, the QNF and quan- 
tal fiuctuations caused by coupling to SCQRPA vibrations are included into the equations 
for the pairing gap and particle number. Under the infiuence of these SCQRPA correc- 
tions, the temperature functional of the quasiparticle occupation number deviates from the 
Fermi-Dirac distribution of independent quasiparticles. The results obtained within the FT- 
BCSH-SCQRPA for the total energies and heat capacities agree fairly well with the exact 
solutions of the Richardson model [1^, [l^ at finite temperature, and those obtained within 
the finite-temperature quantum Monte Carlo method for the realistic ^^Fe nucleus 15 1. 

The positive results of the FTBCSl +SCQRPA encourage a further extension of this 
approach to include the effect of angular momentum on nuclear pairing so that it can be 
applied to study hot rotating nuclei. The rotational phase of nucleus as a whole, such as that 
present in spherical nuclei, or that about the axis of symmetry in deformed nuclei, is known 
to affect nuclear level densities. The relationship between this noncollective rotation and 
pairing correlations has been the subjects of many theoretical studies. The effect of thermal 
pairing on the angular momentum at finite temperature was first examined by Kammuri in 



Ref. 



161 ]. who included in the FTBCS equations the effect caused by the projection M of 



the total angular momentum operator on the 2;-axis of the laboratory system (or nuclear 
symmetry axis in the case of deformed nuclei). It has been pointed out in Ref. IGj that, 
at finite angular momentum, a system can turn into the superconducting phase at some 
intermediate excitation energy (temperature), whereas it remains in the normal phase at low 



and high excitation energies. This effect was later confirmed by Moretto in Refs. [17|, [ISj by 
applying the FTBCS at finite angular momentum to the uniform model. It has been shown 
in these papers that, apart from the region where the pairing gap decreases with increasing 
both temperature T and angular momentum M, and vanishes at a given critical values Tc and 
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Mc, there is a region of M , whose values are shghtly higher than Mc, where the pairing gap 
reappears at T = Ti, increases with T at T > Ti to reach a maximum, then decreases again 
to vanish at T > T2. This effect is called anomalous pairing or thermally assisted pairing 
correlation. In the recent study of the projected gaps for even or odd number of particles 
in ultra-small metallic grains in Ref. [19|] a similar reappearance of pairing correlation at 



finite temperature was also found, which is re: 
this phenomenon was further studied in Refs. 



'erred to as the reentrance effect. Recently, 



20 



2l| by performing the calculations using 



the exact pairing eigenvalues embedded in the canonical ensemble at finite temperature 



and rotational frequency. The results of Refs. [20|, l2l|| also show a manifestation of the 
reentrance of pairing correlation at finite temperature. However, different from the results 
of the FTBCS theory, the reentrance effect shows up in such a way that the pairing gap 
reappears at a given T = Ti and remains finite at T > Ti due to the strong fluctuations of 
the order parameters. 

The aim of the present study is to extend the FTBCSl (FTBCSl+SCQRPA) theory of 
Ref. jll| to finite angular momentum so that both the effects of angular momentum as well 
as QNF on nuclear pairing correlation can be studied simultaneously in a microscopic way. 
The formalism is applied to a doubly degenerate equidistant model with a constant pairing 
interaction G and some realistic nuclei, namely ^''O, ^^Ne, and '^^Ca. 

The paper is organized as follows. The FTBCSl+SCQRPA theory is extended to include 
a specified projection M of the total angular momentum on the axis of quantization in 
Sec. [TTl The results of numerical calculations are discussed in Sec. IIIII The last section 
summarizes the paper, where conclusions are drawn. 



II. FORMALISM 



A. Model Hamiltonian 



We consider a system of particles interacting via a pairing force with the parameter 
G, and rotating about the symmetry axis (noncollective rotation) at an angular velocity 
(rotational frequency) 7 with a fixed projection M (or K) of the total angular momentum 
operator along this axis. For a spherically symmetric system, it is always possible to make 
the laboratory-frame z axis, taken as the axis of quantization, coincide with the body-fixed 
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one, which is ahgned within the quantum mechanical uncertainty with the direction of the 
total angular momentum, so that the latter is completely determined by its 2;-axis projection 
M alone. As for deformed systems, where the axis symmetry is the principal (body-fixed) 
axis, this noncollective motion is known as the "single-particle" rotation, which takes place 
when the angular momenta of individual nucleons are aligned parallel to the symmetry axis, 
resulting in an axially symmetric oblate shape rotating about this axis. Such noncollective 



motion is also possible in high-i^ isomers 



22], which have many single-particle orbitals 



near the Fermi surface with a large and approximately conserved projection K of individual 
nucleonic angular momenta along the symmetry axis. Therefore, without losing generality, 
further derivations are carried out below for the pairing Hamiltonian of a spherical system 



rotating about the z axis 



3, 



ItI . isl - namely 



H = Hp - XN --fM , (1) 
where Hp is the well-known pairing Hamiltonian 

Hp = J2 'k{Nk + N^k) -GY, PlPk' , iV±fe = 4,a±fe , Pi = 4al, , (2) 

h h h' 

with a^^ (a±fe) denoting the operator that creates (annihilates) a particle with angular 
momentum k, spin projection ±mfc, and energy efc. For simplicity, the subscripts k are 
used to label the single-particle states \k, rrik) in the deformed basis with the positive single- 
particle spin projections m^, whereas the subscripts —k denote the time-reversal states 
\k,—mk) {rrik > 0). The particle number operator N and angular momentum M can be 
expressed in terms of a summation over the single-particle levels: 

N = Y,iNk + N^k) , M = Y,^k{Nk-N_k) , (3) 

k k 

whereas the chemical potential A and angular velocity 7 are two Lagrangian multipliers to 
be determined. For deformed and axially symmetric systems, the 2;-projection M and spin 
projections should be identified with the projection K along the body-fixed symmetry 
axis and spin projections respectively, which are good quantum numbers jisl . 

By using the Bogoliubov transformation from particle operators, a\ and a^-, to quasipar- 
ticle ones, at and 



ana a^. 



a\ = Ukal + Vka-k , a_fc = MfcCt-fc - Vkal , (4) 



the Hamiltonian ([T]) is transformed into the quasiparticle Hamiltonian as 

+ hkk'iAlAl + Ak'Ak) + J2 ^kk'MkMi 



kk' 



(5) 



kk' 



kk' 



where A/^" and Af/^ are the quasiparticle-number operators, whereas A^ and At are the 
creation and destruction operators of a pair of time-conjugated quasiparticles, respectively: 



They obey the following commutation relations 

[Ak , Al] = 5kk'Vk , where Vk = l- Uk , 

Wk ' A'] = ^kk'Al, , [U^ , Ak'] = -6kk'Ak' . 

The coefficients in Eq. ([5]) are given as 

bk =bkT jruk = (efc - A)(ufc - vl) + 2GukVkY^Uk'Vk' + Gvl =F 7mfc 



(6) 
(7) 

(8) 
(9) 

(10) 



whereas the expressions for the other coefficients a, bk, Ck, dkk', gk{k'), hkk', and qkk' in Eqs. 
§1 and ([lOD can be found, e.g., in Refs. jjj, M M- 



B. Gap and number equations 

We use the exact commutation relations ([8]) and IQ, and follow the same procedure 
introduced in Ref. llj, which is based on the variational method 

d{n) d{n) dvk 



+ 



m,At]) = o 



(11) 



duk dvk duk 

to minimize the expectation value Ti of the pairing Hamiltonian ([5]) in the grand canonical 
ensemble, 

Tr[C'e-^«l 



{O) 



(12) 



with (p) denoting the ensemble or thermal average of the operator O. The following gap 



equation is obtained, which formally looks like the one derived in Refs. 

G 
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12l |. namely 



Afc = -j:^Y{T>kVk')uk'Vk' 



(13) 



Here 

with the quasiparticle energies Ek defined as 

Ek = y(4^^;^i^A)^+A^ , (15) 
where are the renormahzed single particle energies: 

4 = '^ + ^ E(^'' - ^^') (i^l^Uk) + (Am) ■ (16) 

Notice that the diagonal elements (AIaI) are excluded from all calculations because of 
the Pauli's principle. The Bogoliubov's coefficients, Uk and Vk, in Eq. (fT4|) as well as the 
quasiparticle energy Ek in Eq. ( fTSl) contain the self-energy correction —Gvl- It describes 
the change of the single-particle energy ek as a function of the particle number starting 
from the constant Hartree-Fock single-particle energy as determined for a doubly-closed 
shell nucleus. This self-energy correction is usually discarded in many nuclear structure 
calculations, where experimental values or those obtained within a phenomenological po- 
tential such as the Woods-Saxon one are used for single-particle energies, on the ground 
that such self-energy correction is already taken care of in the experimental or phenomeno- 
logical single-particle spectra. As all calculations in the present paper use the constant 
single-particle levels, determined at T = within the schematic doubly-folded multilevel 
equidistant model and within the Woods-Saxon potentials, we also choose to neglect, for 
simplicity, the self-energy correction —Gvl i^om the right-hand sides of Eqs. (n3|) and (fT5|) 
in the numerical calculations. 

The expectation values {AIa\,) and {AlAk') in Eq. f|T6l) are called the screening fac- 
tors. They are calculated by coupling to the SCQRPA vibrations in the next section. The 
quasiparticle-number fiuctuation (QNF) is included into the gap equation following the exact 
treatment: 

{VkVk') = {Vk){Vk') + SATkk' , with SAfkk' = (A4A4') - (A4)(A4') • (17) 
The term 5Mkk' can be evaluated by using the mean-field contraction as 

5Nkk' ^ 5Ml5kk' , (18) 
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with 



= {6U^f + (5Ar-)^ {5Mtf ^ - n±) , 



(19) 



being the QNF for the nonzero angular momentum. The quasiparticle occupation numbers 
are defined as 



- (Af^) . (20) 

From here, one can rewrite the gap equation (fT3|) as a sum of a level-independent part, A, 
and a level-dependent part, ^A^, namely 



Afe = A + 5Ak 



(21) 



where 



with 



A = GY,Uk'Vk'{Vk') , 6Ak = G 

k' 



UkVk 



(22) 



(V,) = l-nt-n, . (23) 

Within the quasiparticle mean field, the quasiparticles are independent, therefore the 
quasiparticle-occupation numbers (12U]) can be approximated by the Fermi-Dirac distribu- 
tion of non-interacting fermions in the following form 

1 



n 



± 



(24) 



exp[(^{EkTirnk)] + 1 ' 

The equations for particle number and total angular momentum are found by taking the 
average of the quasiparticle representation of Eq. ([3]) in the grand canonical ensemble (fT2l) . 
As the result we obtain 



1 



(P,) + -(1-(P,)) 



M={M) = J2mkint-n^) . 



(25) 



(26) 



We call the set of equations (12T1) . fl25l) and fl26|l as the FTBCSl equations at finite angular 
momentum. By neglecting the QNF (fT9l) . as well as the screening factors {A\aI.,) and 
{AlAk'), i.e. setting = in Eq. ( |T6l) . one recovers from Eqs. ( 12T]) . (l25i) and f[26l) the 



well-known FTBCS equations at finite angular momentum presented in Refs. 



16 



18|. 
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C. Coupling to the SCQRPA vibrations 



1. SCQRPA equations and screening factors 

The derivation of the SCQRPA equations at finite temperature and angular momentum 
is carried out in the same way as that for T = 0, and is formally identical to Eqs. (46), 



(56) and (57) of Ref. 12|]. The only difference is the expressions for the screening factors 



{AlAli) and {A^Ak') at the right-hand side of Eq. f|T6|) . which are now the functions of 
not only the SCQRPA amplitudes, but also of the expectation values {QtQ^') and {QtQti) 



of the SCQRPA operators. As the details of the derivation are given in Ref. 
expressions are quoted below. The screening factors are given as 



111 ], only final 



fj. /Ill' ^ 

... . (44,) -JLy^K: . E (.^(2*4.) ^ -s^s'fi.)) . (28) 



where 

4 = -^, A = [4r, ui:i:: = xi:^!:: +yi:,yi:' , z^,: = xi:yj::+y^'x{:,, {29) 

with Xjf and y^ being the amplitudes of the SCQRPA operators ^ 

Qi = Y.^'^k'A - yi^k) , = [Qi]' . (30) 

k 

The expectation values of (Qj^Q^^') and (QJ^QJ,/) are found as 

{QlQ,') = E y'kyi:' + E(^^T''^..' - i/^^o , (31) 

k kk' 

(QlQl') = - E y'k^i^' + H^u^k^ykk' - wi:i::x,,,) , (32) 

k kk' 

where 

w.^: = x^yi + y^^x^ . (33) 



^ Although at finite angular momentum, the expectation value ([Sfe,S|./]) = Skk'{n^ ~ ^t)' where = 
a^a_fc, becomes finite as 7ij7 ^ n^, the scattering operators Bk and Bj. do not contribute to the QRPA 
because they commute with operators A]., Ak, and 7\4 
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From Eqs. (127|) . (1281) . (131]) and (132|) . the set of exact equations for the screening factors is 
obtained in the form 

= E + E (f^i^^''^^ - , (34) 

kik'j^fifj.' kik[ fMf^' 

= E ^^'-^^^ + E y'" - f/.T>.'^;) • (35) 

H k"fj,fj.' 
2. Quasiparticle occupation numbers 

The quasiparticle occupation numbers ( l20l) are calculated by coupling; to the SCQRPA 



26| . By representing 



phonons making use of the method of double-time Green's functions 
the Hamiltonian (jSj) in the effective form as 



k —k k' pi kfj, 

with Li)fj_ denoting the phonon energies (eigenvalues of the SCQRPA equations) and the vertex 
given as 

= E^^^(^')vW('^.'^ + 3^."') , (37) 

k' 

we introduce the following double-time Green's functions for the quasiparticle propagations 



as well as those corresponding to quasiparticle(8>phonon couplings 



(38) 



T-^,,{t-t') = ((«±fc(t)Q,(t);4,(0)) , T-^tit-t') = ((«±.(t)Qi(t);4,(t'))) • (39) 



Following the same procedure in Ref . HI] , we obtain the final equations for the quasiparticle 
Green's functions G±k{E) in the following form 

1 1 



(40) 
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where 



E 



E^ + u.Y + e' 



Era' 



[UJ 



+ 



Et + ^m)' + 



(41) 
(42) 
(43) 

(44) 



In Eqs. (H2l) - (jSj), the imaginary parts '^^{oj) {uj real) of the analytic continuation of 
M^{E) into the complex energy describe the damping of quasiparticle excitations due to 
coupling to SCQRPA vibrations, = (Q^Q^) is the phonon occupation number, and e is 
a sufficient small parameter. These results allow to find the spectral intensities J^{uj) from 
the relations J^{uj) = i[G±k{uJ + ie) — G±k{uJ — z£:)]/[exp(/5co') + 1] in the form 

1 7.^(^)(e'^- + l)-^ 



7^\u-Et-M^{u)Y + Yit{^)f 
and, finally, the quasiparticle occupation numbers (!20l) as 



Jj^ {u)duj . 



(45) 



(46) 



In the limit of quasiparticle damping 7^(u;) 0, nf can be approximated with the Fermi- 
Dirac distribution 

(47) 



exp(/5E±) + 1 



where E^ are the solutions of the equations for the poles of the quasiparticle Green's func- 



tions G±kiuj) (HOl) . namely 

Ei-E^-Mt{E^)=0. (48) 

The particle-number violation inherent in the BCS-based theories still causes some quan- 
tal fiuctuation of particle number starting from T = 0. This defect can be removed by 
carrying out a proper particle- number projection (PNP). Among different methods of PNP, 
the Lipkin-Nogami (LN) prescription (LN-PNP) 27| is widely used because of its simplicity. 
This method has been implemented into the FTBCSl and FTBCSl+SCQRPA in Ref. (uj, 
and the ensuing approaches are called the FTLNl and FTLNl+SCQRPA, respectively. 
Their extension to M 7^ is straightforward. It is easy to see that, in the nonrotating limit 
(7 = 0), one has = = bk from Eq. (ITO!) . nt = 'f\ from Eqs. P6|) . and all above-derived 



formalism reduces to that presented in Ref. 



ll|. 
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III. NUMERICAL RESULTS 



A. Ingredients of numerical calculations 

The numerical calculations are carried out within a schematic model as well as several 
nuclei with realistic single-particle spectra. For the schematic model, we use the one with 
particles distributed over Q = N doubly-folded equidistant levels. These levels interact via 
an attractive pairing force with the constant parameter G. When the interaction is switched 
off, all the lowest Q/2 levels are filled up with particles so that each of them is occupied 
by two particles with the spin projections equal to im^ (fc=l,. . . , fl, and = 1/2, 3/2, ... 
, Q — 1/2). The single-particle energies are measured from the middle of the spectrum as 
ek = e[k — {Q + l)/2] so that the energies of the lower 5 levels are negative, whereas those 
of the upper 5 levels are positive. The results obtained for = 10, e = 1 MeV, and G = 
0.5 MeV are analyzed in the present paper. 

As for the realistic nuclei, we carry out the calculations for neutrons in and ^'^Ca, 
whereas the contribution of proton and neutron components to nuclear pairing is studied 
for the well-deformed ^^Ne nucleus, where a backbending of moment of inertia as a function 
of the square of angular velocity was detected [28|. The calculations use the single-particle 



energies generated at T = within deformed Woods-Saxon potentials. For the slightly 
axially deformed ^°0, the multipole deformation parameters P2, Ps, Pi, P5, and Pq are 
chosen to be equal to 0.03, 0.0, -0.108, 0.0, and -0.003, respectively. For ^^Ne, the axial 
deformation is rather strong with these parameters taking the values equal to 0.326, 0.0, 
0.225, 0.0, and 0.011, respectively. For the spherical ^"^Ca, all the deformation parameters 
Pi are set to be equal to zero. Other parameters of Woods-Saxon potentials are taken from 



Table 1 of Ref. 29|]. The neutron single-particle spectrum for includes all levels up to 
the shell closure with = 20 (between around -25.84 MeV and 0.49 MeV), from which two 
orbitals, 1(^3/2 and ldi/2, are unbound. These unbound states have been shown to have a 



large contribution to pairing correlations in 20-22Q jgotopes 301]. The neutron single-particle 
spectrum for ^^Ca include all bound states between around -35.6 MeV and -1.05 MeV, up 
to the 2pi/2 orbital of the closed shell with A^ = 50. The single-particle spectra for ^^Ne 
consist of all 11 proton bound states between -30.23 < Cp < -0.156 MeV, and 12 neutron 
ones between -29.834 < e„ < -0.742 MeV. The values of pairing interaction parameter G 
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are chosen so that the pairing gaps A(T = 0, M = 0) obtained at zero temperature and 
zero angular momentum match the experimental values extracted from the odd-even mass 
differences for these nuclei, namely, A(0, 0) ~ 4 MeV for protons in ^^Ne, and 3, 2, and 3 
MeV for neutrons in ^''O, ^^Ca, ^^Ne, respectively. 

The numerical calculations are carried out within the FTBCS and FTBCSl for the level- 
weighted pairing gap A = A^/fi as functions of temperature T, angular momentum M, 
and angular velocity 7. The effect caused by coupling to SCQRPA vibrations is analyzed by 
studying the total energy E = {H) and heat capacity C = dS/dT, whereas the backbending 
is studied by considering the momentum of inertia as a function of 7^ as T varies. 



B. Results within the doubly-folded multilevel equidistant model 

Shown in Figs. [H (a) and[T] (d) are the level- weighted pairing gaps A as functions of T at 
various M, whereas the dependence of A on M at several T is displayed in Figs. [H (b) and 
[T][e). Finally, Figs. [U (c) and[T] (f) show the gaps A as functions of the angular velocity 7 at 
various T. All the results are obtained for the system with = 10 and G = 0.5 MeV, from 
which the left panels are the predictions by the FTBCS theory, whereas the right panels are 
those by the FTBCSl one. It is clearly seen from Figs. [1] (a) and [T] (b) that the FTBCS 
gap decreases with increasing T (M) at M = (T = 0) up to a certain critical value Tc = 
0.77 MeV (Mc = 8h), where the FTBCS gap collapses. The collapse of the pairing gap at 
M = Mc (at T = 0) was proposed by Mottelson and Valatin as being caused by the Coriolis 
force, which breaks the Cooper pairs j3l|. This feature remains with the FTBCS gap as a 
function of M, when T 7^ 0, but with decreasing Mc(T) < Mc as T increases beyond 0.6Tc. 
As for the behavior of the FTBCS g function of T, one notices that, at M slightly 

larger than Mc, the so-called thermally assisted pairing correlation takes place, in which the 
pairing gap is zero at T < Ti, increases at T > Ti to reach a maximum, then decreases 
again to vanish at T > T2 [See. Fig. [1] (a) for M/Mc > 1]. This interesting phenomenon 
was predicted and explained, for the first time, by Moretto in Refs. [l3, Q by applying the 
FTBCS to the uniform model. At M/Mc > 1.1, no FTBCS pairing gap remains. 

Different from the FTBCS predictions, the results obtained within the FTBCSl include 
the effect caused by the QNF. As one can see in Fig. [T] (d), while, in the region of low 
temperature T < Tc, the FTBCSl and FTBCS gaps for different M are rather similar. 
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FIG. 1: (Color online) Level-weighted pairing gaps A as functions of T at various M [(a), (d)], 
and as functions of M [(b), (e)] and 7 [(c), (f)] at several T for = 10, G = 0.5 MeV obtained 
within the FTBCS (left) and FTBCS 1 (right). 

they are qualitatively different at T > T^. Here, the QNF, which is rather strong at high T, 
causes a monotonous decrease of the FTBCSl gap A as T increases. This FTBCSl gap never 
collapses even at very high T. Instead all the values of the FTBCSl gap obtained at various 
M seem to saturate at a value of around 2.25 MeV at T > 5 MeV. This feature shows that, 
the effect of angular momentum on reducing the pairing correlation is significant only at low 
T. In the high temperature region, the QNF leads to a persistence of the pairing correlation 
in a rotating system. Compared to the FTBCS theory, when the QNF is neglected, the 
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effect of thermally assisted pairing correlation also takes place at M/M^ > 1.1. However, 
the FTBCSl gap is now zero at T < Ti, reappears at T = Ti, and remains finite at T > Ti. 
This result is found in qualitative agreement with those obtained in the calculations of the 



canonical gap of in Ref. 



201 ]. where the reappearance of the pairing gap at finite T and 7 



is related to the strong fluctuations of order parameter in the canonical ensemble of small 
systems such as metal clusters and nuclei. In the present paper, we point out the QNF 
as the microscopic origin of this effect. Comparisons between the FTBCSl gaps and the 
canonical ones obtained at (T 7^ 0, M = 0) and (T = 0, M and/or 7 7^ 0) are discussed 
later, in Sec. IIII El 

The QNF has a similar effect on the behavior of the pairing gap A as a function of 
angular momentum. As low T, when the QNF is still negligible, the FTBCS and FTBCSl 
gaps as functions of M are similar. They both decreases as M increases, and collapse at 
M = Mf. and at M slightly higher than Mc for < T /T^ < 0.2, contrary to the trend within 
the FTBCS theory, where Mc(T) decreases as T/Tc increases above 0.6 discussed above 
[Compare Figs. [U (b) and[I] (e)]. At T/Tc = 0.8, e.g., the collapsing points of the FTBCS 
and FTBCSl gaps are M/Mc — 0.85, and 2.9, respectively. 

The FTBCS and FTBCSl pairing gaps are displayed in Figs. [11(c) and[T](f) as functions of 
angular velocity 7 at various T. For T/Tc < 0.2, the pairing gap undergoes a backbending, 
which will be discussed in the Sec. IIII D[ At T/T^ > 0.2 no backbending is seen for 
the pairing gaps. This result agrees with those obtained in calculations within the finite- 
temperature Hartree-Fock-Bogoliubov cranking (FTHFBC) theory, which is applied to the 
two-level model in Ref. 33|]. Within the FTBCSl, the pairing gaps at large M become 
enhanced with T, in agreement with the results obtained within an exactly solvable model 
for a single /7/2 shell in Ref. 21 1. 



C. Results for realistic nuclei 

Shown in Fig. [2] are the level- weighted pairing gaps as functions of T, M and 7 obtained 
within the FTBCS and FTBCSl theories for neutrons in ^'^O. The values of T^ (at M = 0) 
and Mc (at T = 0) are found equal to 1.57 MeV and Ah, respectively. Compared to the case 
with schematic model discussed in the previous section, the difference is that no thermally 
assisted pairing correlation appears within the FTBCS for ^°0. All the FTBCS gaps behave 
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FIG. 2: (Color online) Same as Fig. [T]but for neutrons in ^'^O using G = 1.04 MeV. 

similarly as functions of T with increasing M. At a given value of M, they decrease with 
increasing both T, and collapse at some values Tc(M) < T^. A similar behavior is seen 
for the gaps as functions of M at a given value of T. Here the critical value Mc(T) for the 
angular momentum, at which the gap collapses is found decreasing with increasing T so that 
Mc(T) < Mc [See Figs. [2] (a) and [2] (b)]. Meanwhile, the temperature dependence of the 
FTBCSl gap in Fig. [2](d) shows a clear manifestation of the thermally assisted pairing gap. 
As M increases up to M/M^ — 0.8, the gap decreases monotonously with increasing T up 
to T ^ 1.5Tc, higher than which the gap seems to be rather stable against the variation of 
T. At M/Mc > 0.9, the reentrance of thermal pairing starts to show up as the enhancement 
of the tail at T > Tc. When M/Mc becomes equal to or larger than 1, the gap completely 
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vanishes at low T, but reappears starting from a certain value of T, above which the gap 
increases with T and reach a saturation at high T. At T > 3 MeV, all the gaps obtained at 
different values of M seem to coalesce to limiting value around 0.7 - 0.8 MeV. At a given 
value of T in the region T /T^ < 0.7, as shown in Fig. [2] (e), the pairing gaps decrease steeply 
with increasing T and all collapse at the same value M^. This difference compared to the 
FTBCS theory comes from the presence of the QNF. At T/Tc > 0.8, the QNF becomes 
stronger, which pushes up the collapsing point to Mc(T) = 2Mc. One can also sees some 
oscillation occurring in the region between 0.8 < M/M^ < 1.4 because of the shell structure. 
The collapsing point might be shifted even further to higher M with increasing T, but at 
too high T the temperature dependence of single-particle energies becomes significant so 



3^. 



that the use of the spectrum obtained at T = is no longer valid 

The pairing gaps as functions of angular velocity 7 obtained at various T within the 
FTBCS and FTBCSl theories are plotted in Figs. [2] (c) and [2] (f), respectively. As E^, 
7 and TUk are positive, at T = 0, the quasiparticle occupation number is always zero, 
whereas is a step function of Ek — 'ymk, which is zero if Ek > jmh and 1 if ii^fc < jnik- As 
the result, the FTBCS and FTBCSl gaps decrease with increasing 7 in a stepwise manner 
up to a critical value 7c, where they vanish. At T 7^ 0, the Fermi-Dirac distribution replaces 
the step function, which washes out the stepwise manner in the behavior of the gaps as 
functions of the 7. Here again, once can see that, at T/T^ > 0.8, the QNF is so strong that 
the collapse of the FTBCSl gap is completely smoothed out [Fig. [2] (f)]. 

The level-weighted pairing gaps A for neutrons in ^"^Ca shown in Fig. [3] have a similar 
behavior as as functions of T, M and 7 with the values of Tc and Mc are found to be 1.07 
MeV and 8 h, respectively. The thermally assisted pairing gap appears at M/M^ > 1.0 but 
the high-T tail is much depleted due to a weaker QNF in a heavier system compared to that 
in 200. 

The well deformed nucleus ^^Ne has both neutron and proton open shells, therefore the 
gap and two number equations for protons (p) and neutrons (n) are simultaneously solved 
together with one equation for the total angular momentum M = Mp + M„ to obtain the 
pairing gaps Ap and chemical potential Xp for protons, the corresponding quantities, A„ 
and A„, for neutrons, as well as the angular velocity 7 of the entire nucleus [35!]. The level- 
weighted pairing gaps as functions of T at several M obtained for neutrons and protons 
in ^^Ne are shown in Fig. HI The FTBCS neutron gaps become depleted with increasing 
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FIG. 3: (Color online) Same as Fig. [T]but for neutrons in '^^Ca using G = 0.48 MeV. 

M, and completely disappears at M > 3/i. As a function of T, the FTBCS neutron gaps 
decrease as T increases and collapse at Tc(M), which decreases from Tc(M = 0) ~ 1.7 MeV 
to Tc(M = 3^) ^ 1.1 MeV. The FTBCSl gaps obtained at M < 3 ^ never collapse, but 
gradually decrease with increasing T > Tc{M = 0), and remains a finite value of around 0.4 
MeV at T as high as 4 MeV. At M = 4 ^, whereas there is no FTBCS gap, the thermally 
assisted pairing gap appears within the FTBCSl theory at T > 0, increases with T to reach 
a maximum at T ~ 1.5 MeV, then decreases slowly to reach the same high-T limit of around 
0.4 MeV at T ~ 4 MeV. The situation is the similar for the proton pairing gaps, where the 
effect of thermally assisted pairing correlation takes place at M > 8 h with the rather stable 
values of the gap against T > 3 MeV. 
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FIG. 4: (Color online) Level- weighted pairing gaps as functions of T at various M obtained within 
the FTBCS (left) and FTBCSl (right) for neutrons [(a), (c)], and protons [(b), (d)] in 22Ne using 
Gn = 1.0 MeV and Gp = 1.32 MeV. 

D. Backbending 

For an object that rotates about a fixed symmetry axis, its moment of inertia J is 
found as the total angular momentum M divided by the angular velocity 7, i.e. J = 
M/7. The backbending phenomenon is most easily demonstrated by the behavior of as 
a function of the square 7^. This curve first increases with 7^ up to a certain region of 
7^, where the increase suddenly becomes very steep, and the curve even bends backward. 
This phenomenon is understood as the consequence of the no-crossing rule in the region 



of band crossing 



361] . The SN phase transition has been suggested as one of microscopic 



interpretations of backbending [31 1. 



The values of the moment of inertia obtained at various T within the schematic 
model as well as realistic nuclei, is plotted in Fig. [5l In the schematic model, one can see 
in Figs. [5] (a) and [5] (e) a sharp backbending, which takes place at very low temperatures, 
T/Tc < 0.2. As the QNF is negligible in this temperature region, the predictions by the 
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FIG. 5: (Color online) Moment of inertia as a function of the square 7^ of angular velocity 7 
obtained within the FTBCS (left) and FTBCSl (right) at various T for = 10 [(a), e)], neutrons 
in ^'^O [(b), (f)] and ^^Ca [(c), (g)], and the whole ^^Ne nucleus (including both proton and neutron 
gaps) [(d), (h)]. 

FTBCS and FTBCSl theories are almost identical. As T increases, the moment of inertia 
obtained within the FTBCS changes abruptly to reach the rigid-body value, generating a 
cusp, whereas, under the effect of QNF, the value obtained within the FTBCSl theory 
gradually approaches the rigid-body value in such a way that the cusp is smoothed out. 
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While no signature of backbending is seen in the results obtained in ^"O [Figs. [5] (c) and [5] 
(f)] and ^''Ca [Figs. [5] (d) and[5t^g)], backbending can be seen in ^^Ne [Figs. EKd) and [5] (h)] 



at T < 0.4 MeV in agreement with the experimental data reported in Ref. 



28|. 



E. Corrections due to SCQRPA and particle-number projection 

Shown in Fig. [6] are the level-weighted pairing gaps and moment of inertia, obtained 
within the schematic model with = 10, where predictions offered by several approaches, 
namely the FTBCS, FTBCSl, FTBCSl + SCQRPA, FTLNl, and FTLNl + SCQRPA, are 
collected. In the same figure, the canonical gaps aJ.^^ and aJ.^^ obtained at (T 0, M = 0) 
[Fig. M (al)], (T = 0, M ^ 0) [Fig. M (a2)], and (T = 0, 7 ^ 0) [Fig. M (a3)], are also shown 
(See Appendix |X] for the detailed discussion of the canonical results). 

As seen from Figs. El the effect due to the SCQRPA corrections on the pairing gap 
increases with M. At M/Mc < 0.8 it is rather weak, causing only a slight enhancement of 
the gap at 1.2 < T < 2 MeV as compared with the FTBCSl results [Figs. [6] (al) and [6] 
(bl)]. However, it becomes important at M > 1.2Mc [Figs. [6] (cl), [6] (a2) -[6](c2)], or 7 > 
0.2 MeV/h (at T > 0.8Tc) [Figs. [6] (b3) and [6] (c3)]. In particular, the reappearance of the 
thermal gap at M > l.lMc is significantly enhanced by the SCQRPA corrections [Figs. [6] 
(cl) and[6] (c2)]. For the moment of inertia [Figs. [6] (a4) -[6] (c4)], the SCQRPA corrections 
are important only at low T and 7 < 0.25 MeV/h. At T > 1 MeV, the predictions by all 
the approximations for J' saturate to the rigid-body value. 

As compared to the predictions by the FTBCSl and FTBCSl+SCQRPA, the corrections 
due to LN-PNP are important only at low T and M. As the result, the gap is pushed up to 
be closer to the canonical results at T < and M = [Fig. [6] (al)]. This feature is well- 



known and has been discussed within the present approach at M = in Ref. HI]. At M 7^ 
(7 7^ 0), the effect due to LN-PNP is noticeable in the gaps as functions of M (or 7) only 
at [M < 1.2Mc (7 < 0.2 MeV/h), T < 1 MeV], otherwise the FTLNl (FTLNl+SCQRPA) 
results are hardly distinguishable from the FTBCSl (FTBCSl+SCQRPA) ones [Figs. M 
(b2), [6] (c2), [6] (b3), and [6] (c3)]. Consenquently, for the moment of inertia, the LN-PNP 
corrections to the FTBCSl (FTBCSl+SCQRPA) results are important only at (T < Tc, 
7 < 7c) [Figs. E] (a4) -[6] (c4)]. In particular, the results at T = [Fig. [6] (a4)], where the 
BCSl coincides with the BCS, show that, backbending becomes less pronounced within the 
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FIG. 6: (Color online) Level-weighted pairing gap A and moment of inertia J for = 10 with 
G = 0.5 MeV \e = 0.1 MeV in Eqs. (03]) and dMD]. (al) - (cl): A vs temperature T at different 
angular momenta M. (a2) - (c4): Results obtained at different values of T, namely, (a2) - (c2): 
A vs M; (a3) - (c3): A vs angular velocity 7; (a4) - c4): J" vs 7^. The dotted, thin solid, thick 
solid, thin dash-dotted, thick dash-dotted lines are results obtained within the FTBCS, FTBCSl, 
FTBCSl+SCQRPA, FTLNl, FTLNl+SCQRPA, respectively. The sohd hues with circles and 



boxes in (al) and (a3) correspond to two definitions Ac^'' and Ac^'' of the canonical gaps at T = 0, 
respectively (See Appendix [A]). In (a2) the dashed lines connecting the discrete values of the 
corresponding canonical gaps at T = are drawn to guide the eye. 



.(2) 
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SCQRPA and LN-SCQRPA. For this reason, the corrections due to LN-PNP are omitted in 
the results obtained for reahstic nuclei below. 

Shown in Figs. [7] and [S] are the pairing gaps, total energies and heat capacities as functions 
of T obtained at M/M^ = 0, 0.4, and 0.8 within the FTBCS, FTBCSl and FTBCSl + 
SCQRPA for realistic nuclei, and ^^Ca. The SCQRPA corrections are significant for the 



total energy in the lig^ 
factors (EZD and M. 



it nucleus, ^°0, due to the important contributions of the screening 
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In medium ^"^Ca nucleus, the effect of SCQRPA corrections 
on the total energy is weaker. The corrections due to LN particle-number projection have 
a similar effect as that discussed above for the schematic model, but with much reduced 
magnitudes, so they are not shown in these figures. With increasing M the pairing gap 
decreases. As the result, the total energy becomes larger but the relative effect of the 
SCQRPA correction does not change. For the heat capacity, as has been reported in Ref. 



ll| . the spike at obtained within the FTBCS theory, which serves as the signature of the 



sharp SN phase transition, is smeared out within the FTBCSl theory into a bump in the 
temperature region around Tc. With increasing M, this bump becomes depleted further. 
Finally, the SQRPA corrections erase all the traces of the sharp SN phase transition in the 
model CclSG clS well as realistic nuclei. 



IV. CONCLUSIONS 

The present work extends the FTBCSl (FTBCSl + SCQRPA) theory to finite angu- 
lar momentum to study the pairing properties of hot nuclei, which rotate noncollectively 
about the symmetry axis. The FTBCSl theory includes the quasiparticle number fluctua- 
tion whereas the FTBCSl + SCQRPA also takes into account the correction due to dynamic 
coupling to SCQRPA vibrations. The proposed extension is tested within the doubly de- 
generate equidistant model with = 10 particles as well as some realistic (spherical and 
deformed) nuclei, ^°0, ^^Ne, and ^^Ca. The numerical calculations were carried out within 
the FTBCS, FTBCSl, and FTBCSl + SCQRPA for the pairing gap, total energy, and heat 
capacity as functions of temperature T, total angular momentum M, and angular velocity 
7. The corrections due to the Lipkin-Nogami particle-number projection are also discussed. 
The analysis of the numerical results- allows us to draw the following conclusions: 

1. The proposed approach is able to reproduce the effect of thermally assisted pairing 
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FIG. 7: (Color online) Level-weighted pairing gaps A, total energies i5, and heat capacities C as 
functions of temperature T for three values of angular momentum M obtained within the FTBCS 
(dotted hues), FTBCSl (thin sohd lines) and FTBCSl + SCQRPA (thick sohd hues) for neutrons 
in with G = 1.04 MeV (e = 0.1 MeV). 

correlation that takes place in the schematic model within the FTBCS theory, according to 
which a finite pairing gap can reappear within a given temperature interval, Ti < T < T2 
(Ti > 0), while it is zero beyond this interval. However, this phenomenon does not show up 
in realistic nuclei under consideration. 

2. Under the effect of QNF, the paring gaps obtained within the FTBCSl at different 
values M of angular momentum decrease monotonously as T increases, and do not collapse 
even at hight T in the schematic model as well as realistic nuclei. The effect of thermally 
assisted pairing correlation is seen in all the cases, but in such a way that the pairing gap 
reappears at a given Ti > and remains finite at T > Ti, in qualitative agreement with the 



canonical results of Ref. 



2Q|. 



3. The backbending of the moment of inertia is found in the schematic model and in 
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FIG. 8: (Color online) Same as in Fig. [7] but for neutrons in ^^Ca with G = 0.48 MeV (e = 0.1 
MeV). 

^^Ne in the low temperature region, whereas it is washed out with increasing temperature. 
This effect does not occur in and '^^Ca, in consistent with existing experimental data 
and results of other theoretical approaches. 

4. The effect caused by the corrections due to the dynamic coupling to SCQRPA vibra- 
tions on the pairing gaps, total energies, and heat capacities is found to be significant in the 
region around the critical temperature of the SN phase transition and/or at large angular 
momentum M (or angular velocity 7). It is larger in lighter systems. As the result, all 
the signatures of the sharp SN phase transition are smoothed out in both schematic model 
and realistic nuclei. The SCQRPA corrections also significantly enhance the reappearance 
of the thermal gap at finite angular momentum. On the other hand, the effect caused by 
the corrections due to PNP is important only at temperatures below Tc, and at quite low 
angular momentum. In particular, it makes backbending less pronoucned at T = 0. 

Still the fluctuations due to violation of angular-momentum conservation are not imple- 
mented in the present extension. We hope that further studies in this direction will shed 
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FIG. 9: (Color online) (a) Canonical moment of inertia vs 7^; (b): Absolute values \{£)c ~ £m.i.\ 
(solid line) and |func. I (dotted line) vs 7; (c): [a[1^]^ vs 7 for the schematic model with = 10 
and G = 0.5 MeV at T = 0. 

light on this issue. 
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APPENDIX A: ON THE COMPARISON WITH CANONICAL RESULTS 



It has been shown in Sec. [TTlBlthat the FTBCSl (FTBCSl+SCQRPA) produces results 
in qualitative agreement with the canonical ones of Ref. 201], in particular, the reappearance 
of the thermal gap at M 7^ 0. However, it is important to make clear the difference between 
the predictions by BCS-based approaches and the canonical results. As a matter of fact, 
the z-projection M of the total angular momentum within the FTBCS (FTBCSl) approach 
is temperature-independent. At T varies, by solving the FTBCS (FTBCSl) equations, the 
angular velocity 7 is defined as a Lagrangian multiplier so that M, being the thermal av- 
erage of the total angular momentum within the grand canonical ensemble (JT2l) . remains 
unchanged. In this way, within the FTBCS (FTBCSl), the angular velocity 7 varies with 
T, whereas M does not. Similar to that for choosing the chemical potential A to preserve 
the (grand-canonical ensemble) average particle-number N, this contraint is physically rea- 
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sonable when the total angular momentum is conserved as in the noncollective rotation of 
spherical systems or rotation of axially symmetric systems about the symmetry axis, as has 
been discussed in Sec. Ill A[ 

On the contrary, the canonical results in Ref. 20|] are obtained by embedding the eigen- 



values Ey^i{'~f) = Ey — '-fMi.^i in the canonical ensemble with the partition function 

Z(/3,7) = ^e-^^-W . (Al) 

Here Ejy denote the eigenvalues of the z/th state with seniority z/ at 7 = 0, whereas M^^i 
are the z-projections of angular momenta of u nucleons. While the eigenvalues E^, are 
obtained by separately diagonalizing the pairing Hamiltonian Hp in Eq. ([2]), the rotational 
)art = '^■exp^P'jM^^i) of the partition function Z(/5,7) is calculated following Ref. 



32]. The resulting canonical average value (M(/?,7))c = 7)~^5Z(/3, 7)/97 of angular 

momentum, therefore, varies with T. On the other hand, the angular velocity 7 just plays 
the role of an independent parameter, therefore, does not depend on T. By the same reason, 
each canonical average value {M{j3, 7))c corresponds to a single value of 7, i.e. the canonical 
moment of inertia J^c undergoes no backbending, as shown in Fig. [9] (a). 

Because of this principal difference, a quantitative comparison between the FTBCS (FT- 
BCSl) results, and the canonical ones as functions of M (or 7) at T 7^ unfortunately 
turns out to be impossible. To establish a meaningful correspondence, one needs to know 
the exact eigenvalues of the ground state as well as all excited states of the pairing problem 
described by Hamiltonian ([1]) so that, by embedding the eigenvalues in the grand canonical 
ensemble, 7 becomes a function of T in such a way to keep (M(/3, 7))c always equal to M. 
To our knowledge, this problem still remains unsolved. One may also try to estimate the 
results within the microcanonical ensemble. However, here one faces a problem of extracting 
the nuclear temperature, which is rather ambiguous at low level density (small A^) within 



the schematic model under consideration 
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391], whereas the extension of exact solution 



of the pairing problem to T 7^ is unpracticable at A^ > 16. 

Therefore, in the present paper, we can only compare the predictions of our approach 
with the canonical results as functions of temperature T at M = 0, or as functions of M (or 
angular velocity 7) at T = 0. For this purpose, and given several definitions of the "effective" 
gaps existing in literature, we choose to employ in the present paper two definitions of the 
canonical gaps, A^.^-* and They should be understood as effective ones since a gap 
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per se, which is a mean-field concept, does not exist in the exact solutions of the pairing 
problem. 

The canonical gap A^.^-* is defined from the pairing energy Sp^h of the system as 

l^C 1 ~ ^^pair ) ^pair = (^)c ^^m.f. ^unc. ; ^m.f. = 2 ^ ^ Cfc/fc ; ^unc. = Cr ^ ^ . 

k k 

(A2) 

Here {S)c is the total energy within the canonical ensemble with the partition function 
Z{j3,'-f) given by Eq. flAip of a system rotating at angular velocity 7, or with the partition 
function Z{P, 0) at M = 0. The term S^s. denotes the energy of the single-particle motion 
described by the first term at the right-hand side of the pairing Hamiltonian Hp in Eq. . 
Functions fk are occupation numbers of kth orbitals within the canonical ensemble. The 
energy S^,i, becomes that of the mean-field once the single-particle occupation numbers fk 
are replaced with those describing the Fermi-Dirac distributions of independent particles. 
The energy S^nc. comes from the uncorrelated single-particle configurations caused by the 
pairing interaction in Hamiltonian ([2]). Therefore, by subtracting the term S^j, + Sunc. from 
the total energy {S)c, one obtains the result that corresponds to the energy due to pure 

pairing correlations. The definition ( 1A2I) is very similar to that given in Ref. js^]- It is, 

(2) 

however, different from the canonical gap , which is used in Refs. The latter is 

defined as 

[Ag)]2 = -G[(f)c-(^(G = 0))c] , (A3) 

where {S{G = 0))c is the total canonical energy {£)c at G = 0. 

The canonical gaps A^c^ and Ag^ are shown in Figs. [6](al),[6](a2), and[6](a3) as functions 
of temperature T (at M = 0), angular momentum M (at T = 0), and angular velocity 7 
(at T = 0), respectively. It is seen from these figures that the difference between the two 
canonical gaps A[^^^ and A^^ is rather significant at large T for M = 0, and at large M 
(or 7) for T = 0. The reason is rather simple since the definition flA2p of A[!'' is rather 
similar to that for the BCS gap. As a matter of fact, by replacing the canonical single- 
particle occupation numbers fk with the Bogoliubov's coefficients vl, and the total energy 
{S)c with that obtained within the BCS theory, the gap Aq^ reduces to the usual BCS gap. 
Meanwhile, by doing so with A^-*, the energy {S{G = 0))c just reduces to the Hartree- 
Fock energy, leaving the uncorrelated energy £^unc. out of the definition. Consequently, as 
functions of T, the gaps predicted by the BCS-based approaches under consideration agree 
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better with the canonical gap A[!^ than with A[?-* [Fig. [6] (al)]. 

As functions of angular velocity 7, both the squared values (lA2p and (]A3I) of the canonical 
gaps undergo a stepwise decrease with increasing 7. The step occurs whenever the state of 
the lowest energy changes from u — 2 to u, causing a stepwise increase of (M(/5, 7))c |20 |. 
Therefore, for = 10, the pairs are gradually broken in 5 steps with a corresponding 
stepwise increase of seniority u from to 10 by two units in each step. However, Fig. [9] (b) 
shows that the absolute value of the uncorrelated energy f^uno which enters in the definition 
(IA2I) of the gap A^q\ becomes larger than that of the difference Sq — i^m.f. already at the 

second step, leading to [A^"*]^ < [Fig. [9](c)], i.e. an imaginary value for Aq\ As the result, 

(2) 

instead of collapsing as A^ in 5 steps at a rather large value of M (or 7), the canonical gap 
Aq^ collapses in two steps at a value of M (or 7) much closer to Mc (or 7c) for the BCS 
gap [Figs. M (a2) and [6] (a3)]. Once again, this makes the gaps predicted by the BCS-based 
approaches as functions of M (or 7) agree better with the canonical gap A[^^\ rather than 
with [Figs. [6](a2) and[6](a3)]. 
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